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Abstract 

A Dual Parton Model with a formation zone intranuclear cascade in 
the spectators of the projectile and target nuclei is studied. The hadrons 
produced in the formation zone cascade contribute to Feynman-x^ and 
lab-x distributions in the fragmentation regions of the target and projectile 
nuclei. We discuss the consequences of this model in the secondary cosmic 
ray production, by analyzing the calculated spectrum weighted moments 
for pion and kaon production. We show that the proposed model leads to 
significant differences with respect to a simple superposition model, where 
the nucleus-nucleus collision is replaced by a few corresponding nucleon- 
nucleus collisions. 
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1 Introduction 



Nucleus-nucleus collisions are of great importance for the understanding of the 
cosmic ray cascades in the atmosphere. Models for sampling hadron production 
events in nucleus-nucleus and hadron-nucleus collisions are needed for the simu- 
lation of the development of cosmic ray showers in atmosphere. In some models 
used for the simulation of cosmic ray cascades, examples are HEMAS [|TJ and 
SIBYLL ||, the correct treatment of nucleus-nucleus collisions is replaced by a 
simpler superposition model [], where the nucleus-nucleus collision is replaced by 
the corresponding nucleon-nucleus collisions. Here we study a model where the 
nucleus-nucleus collisions are treated in a more detailed way. 

Soft and hard multiple interactions between nucleons of both nuclei dominate 
the hadron production in most of the kinematic region covered by the interaction 
and are well described in the framework of the Dual Parton Model (DPM) ||. 
They were extensively studied with the Monte Carlo (MC) implementation of 
this model: dpmjet-ii |4]]. This event generator was already applied for sampling 
cosmic ray cascades @, |5], §]. However, when dealing with particle or fragment 
production in the forward or backward fragmentation regions a detailed descrip- 
tion of intranuclear cascade processes and of nuclear disintegration is important. 
A formation zone intranuclear cascade in the target nucleus was therefore consid- 
ered since a few years J7|, || || |J. A formation zone intranuclear cascade (FZIC) 
model in both the target as well as the projectile nucleus, the calculation of nu- 
clear excitation energies, models for nuclear evaporation, high energy fission and 
break-up of light nuclei were discussed for hadron-nucleus and nucleus-nucleus 
collisions in [pi], [11). It was shown that these MC implementations for hadron- 



nucleus and nucleus-nucleus interactions describe successfully the basic features 
of target and projectile associated particle production. 

The model was compared in [10j [l]J mostly to data for asymmetric collisions 
in which projectile and target masses are different. It seems that the performance 
of the model in these collisions is rather good. Asymmetric collisions are the most 
common ones also in cosmic ray cascades. 

Here we will study this model with then aim to find differences with respect 
to the superposition models mentioned above. 

In Section 2 we summarize briefly the main steps of sampling hadron-nucleus 
and nucleus-nucleus interactions within the event generator dpmjet-ii. Fur- 
thermore, we summarize the basic ideas of the FZIC model. In Section 3 we 
calculate Feynman-x^ distributions of pions produced in hadron-nucleus and 
nucleus-nucleus collisions. In Section 4 we study spectrum-weighted moments 
and energy fractions carried away by some kinds of secondaries and we compare 
them with the expectations from superposition models. It is pointed out that the 



1 it must be noticed that SIBYLL takes however into account the correct differences in the 
interaction heights as coming from the nucleus nucleus cross sections. 
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DPM differs from the superposition model. Finally, in Section 5 we summarize 
our results. 



2 The two-component DPM for hadron-nucleus 
and nucleus-nucleus collisions 

2.1 The DPMJET-II event generator 

The two-component DPM and its MC realizations for hadron-hadron, hadron- 
nucleus and nucleus-nucleus collisions have been discussed in detail in [H| fl. 
Therefore, we briefly summarize the main steps leading to the multiparticle state, 
which is the starting point for the intranuclear cascade. 

The MC model for hadron-nucleus and nucleus-nucleus interactions starts 
from an impulse approximation for the nucleons of the interacting nuclei. The 
spatial initial configuration, i.e. the positions of the nucleons in space-time in 
the rest system of the corresponding nucleus, is sampled from standard density 
distributions. For energies above 3-5 GeV/nucleon the collision proceeds via v 
elementary interactions between v v and v t nucleons from the projectile and target, 
respectively. The values u, u p , and v t are sampled according to Glauber's multiple 
scattering formalism using the MC algorithm of ]14|]. The particle production is 



well described by the two-component DPM which is applied as in hadron-hadron 
interactions [[12], [13|]. As a result, a system of chains connecting partons of the 
nucleons involved in the scattering process is formed. The chains are hadronized 



applying the JETSET model [JT5], [L6|]. The hadrons may then cause intranuclear 
cascade processes, which are treated by the FZIC model [fjj, an extension of the 
intranuclear cascade model of ref. |L7|, At energies below 3-5 GeV/nucleon the 



FZIC model is able by itself to describe reasonably the inelastic nuclear collisions. 

2.2 The formation zone cascade in target and projectile 
nuclei 

In the following we summarize the main ideas of the FZIC model for hadron- 
nucleus and nucleus-nucleus interactions 0, ||, ^ [|, (E(J [T| . The physical picture 



explaining the absence of the intranuclear cascade at high energies is the concept 
of the formation zone |L9| . It has been introduced in analogy to the Landau- 
Pomeranchuk [[21J effect, which explains the observation that electrons passing 
through high density materials become more penetrating at high energies. For the 
formation zone of an electron with 4-momentum p and energy E upon radiation 
of a photon with 4-momentum k one obtains 

E E1 a) 



k ■ p muj, 
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where u e is the frequency of the photon in the rest frame of the electron and 
E /mis the time dilatation factor from the electron rest frame to the laboratory. 
Within the quark model, the states being formed in the primary nucleon-nucleon 
interaction can be understood as consisting of valence quarks only (i.e without 
the full system of sea quarks, antiquarks and gluons) and have therefore a re- 
duced probability for hadronic interactions inside the nucleus [[F]. The formation 
zone concept can be translated to hadron production as follows f|. We consider 
the formation zone cascade in the rest system of the target nucleus (laboratory 
system) or in the rest system of the projectile nucleus. Denoting the 4-momenta 
of the projectile hadron p p and of the secondary hadron p s with 



Pp = (E p , 0, 0, yjE% - ml), Ps = (E s ,p sL , sjE* - ml - f sA _) (2) 

and replacing in Eq. (|]) the electron momentum by p p and the photon momentum 
by p s , the hadron formation time is, for E p 3> m p : 



2E S E s . . 

T~Lab = 7 T5 ; 5 ; o ' X = "FT" \ ) 

{m p xy + m z s + p s± E p 



The term (m p x) 2 can be neglected for most of the produced secondaries, so one 
can approximate 

^Lab ~ lsT s , Is = — ■ (4) 

m s 

We define an average formation time to create a complete hadronic state r s in 
the rest system of the secondary hadron s 0, |H : 

T « = T 2 i 2 • 5 

m\ + pj ± 

where r is a free parameter, which has to be determined by comparing particle 
production within the model to experimental data. Typical values are in the 
range from 1 fm/c to 10 fm/c. Here we use r = 4.5 fm/c. This value differs 
somewhat from the values reported in ref. |T^, |TI[. The reason for this is an 
updated treatment of nuclear evaporation with respect to the version of the code 
used in |10| , |IT| . For each secondary we sample a formation time r from an 



exponential distribution [21] with an average value as given in Eq. (y). As it was 
described in 0, in the MC model the full space-time history of the collision is 
known. 

After having assigned a formation time to a secondary, its spatial coordinates 
in the rest system of both nuclei are known and we start with considering an 
intranuclear cascade step in one (randomly chosen) of the spectators. 

Due to relativistic time dilatation, those secondaries which have sufficiently 
high energies in the rest frame of the considered nucleus are formed mostly out- 
side of the spectator part of this nucleus, whereas those with lower energies are 
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formed inside. The latter may penetrate the spectator and initiate intranuclear 
cascade processes. Elastic and inelastic interactions with spectator nucleons are 
treated using the MC- model hadrin ||22|| . This code is based on measured cross 
sections and interaction channels up to a laboratory momentum of 5 GeV. We 
apply hadrin to hadron-nucleon interactions up to 9 GeV and neglect those at 
higher energies. Reinteractions beyond 5 GeV occur much less frequently than 
reinteractions below 5 GeV and a more detailed treatment would not change the 
results discussed in this paper. Furthermore, we take into account absorption of 
low-energy mesons and antiprotons by interactions with two-nucleon systems (for 
pion absorption we use the cross sections as given by Ritchie p3| ) and Pauli's 



principle ||. In case no interaction is possible in the considered spectator, we 
proceed with sampling a cascade step in another spectator. 

For secondaries produced in intranuclear cascade processes we apply the same 
formalism, i.e. a formation time is sampled, the secondary is transported to the 
end of the formation zone and reinteractions are treated if they are possible. 
Due to these intranuclear cascade processes, nucleons are knocked out of the 
residual spectator nuclei if their energy is high enough to escape from the nuclear 
potential. 



3 Feynman-Xf distributions in hadron— nucleus 
and nucleus— nucleus collisions 

The so-called "cascade particles" are the particles which are produced or knocked- 
out off the spectator nucleus by the FZIC. Target associated cascade particles 
are experimentally studied (mostly in emulsion experiments) and are often called 
"grey" prongs. 



In the papers |10], [Ijj most of the available data on grey particle production in 
high energy hadron-nucleus and nucleus-nucleus collisions are compared to the 
model we use. Reasonable agreement is found as far as the average multiplicities 
of grey prongs, their multiplicity and angular distribution, and their correlations 
with the number of fast produced particles are concerned. There are two kinds 
of particles contributing to the grey prongs: protons and charged pions. The 
emulsion experiments are not able to differentiate between such protons and 
charged pions. There is however one heavy ion experiment (WA80 |24], which 
presents only protons in the energy region of the grey prongs in their multiplicity 
distributions. We can only draw some rather indirect support for the ability of our 
model to describe the pions from the intranuclear cascade. This support comes 
from the fact that the model agrees as well to the emulsion data on grey prongs 
as to the WA80 data on protons alone. Roughly 20% of the grey prongs (defined 
with Lorentz parameter 0.23 < (5 < 0.70) in the model are charged pions|jTT] . 



The Feynman-rr^ distributions are the most suitable way to present the parti- 
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cle production in the target and projectile fragmentation regions. Unfortunately, 
most experiments, which measure grey prongs are not able to measure the particle 
momenta and identity and are therefore not able to present Feynman-Xi? distri- 
butions. We define Feynman-Xi? in h-A and A- A collisions with the longitudinal 
momentum in the hadron-nucleon (or nucleon-nucleon) cms p^ ms : 

„cms 

" = j£rr (•> 



-.cms 
z,max I 



where we use a p c ^ ax calculated disregarding Fermi momenta. The Feynman-Xi? 
distribution f(x F ) is defined as follows 

dN 

f{x F ) = x F - — . (7) 

CLXf 

For cosmic ray calculations we are mainly interested in the contributions of 
pions from the intranuclear cascade to the Feynman-Xi? distributions in the frag- 
mentation region of the projectile nucleus, since particles from target fragmen- 
tation are very slow in the lab-frame. However, experimental data about the 
target fragmentation are very important for the tuning of the interaction model. 
The changes in Feynman-Xi? distributions from p-p to p-nucleus collisions in 
the proton fragmentation region have been measured (and have been compared 
to DPMJET-II in 01]), but we are not aware of any such measurements in the 
fragmentation region of the nucleus. Therefore, we have to rely on the model, but 
we stress that it would be highly desirable to measure Feynman-Xi? distributions 
in the fragmentation region of target or projectile nuclei. 

In Fig. [I] we present the Feynman-Xi? distribution of ir + mesons in p-air 
collisions at a lab-energy of 500 TeV . There are two plots, one for the full model 
with formation zone intranuclear cascade in the target nucleus, and one without 
this cascade. We also compare the distributions with the one in p-p collisions 
at the same energy. In the target fragmentation region, at x F values between 
-0.3 and -1, we find significant differences in the distributions due to the pions 
produced by the formation zone cascade in the target nucleus. The difference 
between the two plots is not due to a large number of particles. If we define 
"grey" the particles with Lorentz-/? of 0.23 < (3 < 0.70 in the lab-frame, we 
find in the calculation 0.06 grey 7r + per interaction, on average. However, not all 
charged pions from the formation zone cascade belong to the grey particles. 

Next, we turn to nucleus-air collisions and calculate again Feynman-Xi? distri- 
butions. It is already enough to look at Fig.|l] to understand, qualitatively, what 
happens in nucleus- nucleus collisions. Fig.|I] gives the Feynman x F distribution 
in N-p collisions if we simply exchange x F by —x F . Therefore, we expect that 
the FZIC will lead to significant changes of the distributions in the fragmentation 
regions of nuclear projectiles in nucleus-proton and nucleus-nucleus collisions. 

In Figs. |] to | we present the Feynman-x^ distribution of 7r + mesons at 
a lab-energy of 500 TeV per projectile nucleon in He-air, O-air, and Fe-air 
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collisions. The model, as shown in 0], exhibits a rather good Feynman-scaling 
behaviour, so that these distributions look rather similar also at other different 
energies. Fig.^ is included mainly to demonstrate that for a light nucleus, like 
He, the formation zone cascade makes really no difference and can be neglected. 
There are two plots on each figure, one for the full model with formation zone 
intranuclear cascade in the target as well as in the projectile nucleus, and one 
without this cascade. We also compare the distributions with the corresponding 
p-air collisions at the same energy and with the formation zone cascade. In the 
target fragmentation region, at xp values between -0.3 and -1 we find in each 
plot nearly the same differences in the distributions due to the pions produced 
by the formation zone cascade in the target nucleus. Such differences are also 
visible in the projectile fragmentation region at positive Xp values, but here those 
differences, which were not significant in the case of He-air collisions, are found 
to rise with the mass number of the projectile nucleus. Clearly, the formation 
zone cascade in the spectators of the Fe nucleus produces more pions than the 
one in the O-spectators. 

4 Spectrum weighted moments, energy frac- 
tions and comparison with superposition mod- 
els 



Following for instance the basic discussion of Ref. [2B|, we introduce a variable 



Xi a b similarly to Feynman-Xi?, but this time in the lab-frame : 

Xlab = |k (8) 

Ei is the lab-energy of a secondary particle i and E is the lab-energy of 
the projectile in a h-A collision (or the energy per projectile nucleon in a A-A 
collision). We introduce xi a b distributions F{xi a j ) ) : 

Fi{X lab) = X iab - • (9) 

We note that the Feynman-Xi? distribution at positive xp in the projectile 
fragmentation region is a very good approximation to the xi a b distribution. There- 
fore, the Feynman xp distributions given in the last Section give also a rather 
good picture for xi a b distributions. 

The cosmic ray spectrum weighted moments in A-B collisions are defined as 
moments of the F(xi a b) 

Zf~ B = / {XlabV^Ft^Xia^dXiab- (10) 
J 
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Here —7 ~ -1.7 is the power of the integral cosmic ray energy spectrum. 

The spectrum weighted moments for nucleon-air collisions, as discussed in 



Ref. [26], determine the uncorrelated fluxes of energetic particles in the atmo- 
sphere. This result has been obtained for cosmic ray cascades initiated by a 
hadron primary. There, the spectrum weighted moments (which change only 
slightly with the primary energy) are the same for all generations of the cascade. 
We can not apply this result directly for our situation. A cosmic ray cascade 
initialized by primary nuclei becomes after the first two to three generations also 
a cascade with only hadrons participating. There are no theoretical arguments 
for the relevance of the spectrum weighted moments of nucleus-nucleus collisions 
for the cosmic ray cascade as a whole. These moments should be relevant only 
for the first generations. However, we might use the spectrum weighted moments 
and energy fractions just as any other moment of the Feynman-x^ distribution 
with the aim to point out differences between models. In particular, we use these 
moments to stress that the superposition model for nucleus-nucleus collisions is in 
reality only a modest approximation to the correct treatment of nucleus-nucleus 
collisions. 

We also introduce the energy fraction Kf~ B . In h-A collisions that is the 
fraction of primary energy carried by secondaries of type i; in A-B collisions it is 
this energy fraction multiplied with the mass number A of the projectile nucleus 

K?~ B = I** F A ~ B (xi ab )dx lab , (11) 



n 



where Aa is the mass number of the projectile nucleus A. After recalling the 
concept of superposition model, in the following we shall compare the calculation 
of spectrum weighted moments with our model in different conditions, pointing 
out the differences with respect to a simple superposition model. 



4.1 The superposition model 

Generally, we call superposition model the approximation in which in the cosmic 
ray cascade the collision of a nucleus A, with total energy E, against a target B, is 
treated as the superposition of Aa independent nucleon-B collisions, each nucleon 
having an energy E/ Aa- This model is based on the hypothesis that, when the 
energy per nucleon of the projectile is much larger than the single nucleon binding 
energy, the Aa nucleons will interact incoherently. In the context of cosmic ray 
simulations, for a primary of total energy E and mass number Aa (A^>1), the 
cascade generated is equivalent to the total effect of Aa showers initiated by Aa 
independent nucleons of energy E/ Aa- For example, the cascade for a primary 
iron nucleus of 560 TeV would be simulated as the superposition of 26 proton 
+ 30 neutron initiated showers, each of 10 TeV. Examples of application of the 
superposition model to cosmic ray calculations can be found in Ref. |27|, |28| . 



S 



Here, in the comparison with DPMJET-II we will apply a less general con- 
cept of superposition model which we call restricted superposition model, making 
reference to a single nucleus-nucleus interaction. In a minimum bias A-B colli- 
sion, not all Aa projectile nucleons interact. The number of interacting projectile 
nucleons is only v^~ B '. For the collisions considered the average values of v^~ B 
are always given in Tables 1 to 4. Obviously, in order to make a comparison with 
DPMJET-II, we should use a restricted superposition model where only v^~ B 
projectile nucleons interact. In this restricted superposition model we have for 
instance Z° upA ~ B = u^' B ZP- B (for the same energy /nucleon of the projectile). 



4.2 The Glauber model in nucleus— nucleus collisions and 
the restricted superposition model 

The Glauber model for nucleus-nucleus collisions f29| is used by most models for 



nucleus-nucleus and nucleon-nucleus collisions. DPMJET uses the implementa- 



tion of ref. [p^4[ . This model has properties different from those of the restricted 
superposition model. In order to understand that, we recall the main formulae 
for inelastic nucleus-nucleus collisions. From these expressions we understand 
that the Glauber model cannot be reduced just to nuclear geometry. 

We start with the scattering amplitude of two nuclei with mass numbers A 



and B in the impact parameter representation |3C| gl], ^1 

F(b) A _ B =< V>Ml - II fit 1 - X(b - *j + r k )}\^\ > (12) 

3=1 k=l 

where b is the impact parameter vector and Sj and are the coordinates of 
the nucleons with respect to the centers of the nuclei A and B , respectively, 
in the impact parameter plane. The ip l A B and ip A B are the initial and final 
state wave functions of nuclei A and B, while x(b) is the dimensionless elastic 
nucleon-nucleon scattering amplitude in the impact parameter representation. 
The (Gaussian) parametrization of x(b) used in DPMJET-II was described in 
detail in [Q . 

Furthermore, we use a simple assumption for the squares of the state wave 
functions 

A B 
\^A? = YlpAiSjiZj), \tfj B \ 2 = Y[pB(T k ^k) (13) 

j=l k=l 

where p A and p B are the one particle densities of the nuclei A and B. 
From this we get the total A-B cross section 

af_ B = 2Re J d 2 bF(b) A -B- (14) 

The inelastic A-B cross section is determined by the profile function T(b) 

oTb = I d 2 hT{h) (15) 
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T(b) is defined as follows 



r(b) = /{i-nn(i -p^HUpA^dd^Hll pB^dh,} (16) 

i=l j=l i=l j=l 

where 

Pa = x( b - s * + t j) + x*( b - s i + t j) ~ x( b - s i + r i)x*( b - ^ + rj). (17) 

This profile function, which contains more than the nuclear geometry, is also used 
to sample in the Monte Carlo calculation the impact parameters of the inelastic 
A-B collisions. 

Sampling the impact parameters of nucleon-nucleus and nucleus-nucleus col- 
lisions simply from geometry, will certainly lead to a restricted superposition 
model. Unfortunately, it is somewhat cumbersome to see, without the actual nu- 
merical calculations, that the profile functions according to the Glauber model for 
nucleon-nucleus collisions differ sufficiently from the profile functions for nucleus- 
nucleus collisions to prevent the restricted superposition model to give identical 
results as the full model. As an example related to the case discussed in the next 
subsection, we show in Fig.|5| the profile functions for p-Nitrogen He-Nitrogen and 
Fe-Nitrogen as function of b = |b|. For p-Nitrogen collisions, the profile function 
never becomes completely absorptive (r(6) = 1), even in the center (6=0), while 
in Fe-Nitrogen collisions T(b) =1 is reached for b < 5 fm. 

4.3 Calculation results, DPMJET-II without FZIC, eva- 
poration and determination of the residual nuclei 

In Table |l] and H we present spectrum weighted moments and energy fractions 
as calculated in DPMJET-II in the model without FZIC, evaporation and con- 
sequent determination of the residual nuclei. We consider p-air, He-air, O-air 
and Fe-air collisions at laboratory energies per nucleon of 50 GeV, 500 TeV, and 
5 • 10 6 TeV, the latter energy (5 • 10 18 eV) is the absolute highest energy at which 
DPMJET-II is valid, at present. We also give the average number of projectile 
nucleons u p participating in minimum bias collisions. 

In Tables [l] and |2] we present a further comparison, at the two upper energies, 
of the spectrum weighted moments of pions and kaons and the energy fractions 
of pions, kaons and baryons minus antibaryons according to the model, with the 
moments and energy fractions obtained from the restricted superposition model. 
The baryons are here only the secondary baryons from the collision. Spectator 
baryons from the projectile or target nucleus do not contribute to the energy 
fraction Kr^ B . Since the energy fraction of the newly produced baryons (pro- 
duced in bb pairs) is equal to the energy fraction of the produced antibaryons, 
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the energy fraction K A ~^ B gives the energy fraction retained by the baryons par- 
ticipating in the collision. In the restricted superposition model the spectrum 
weighted moments and energy fractions are given by the simple expressions 

z su P A-B = V A~B Z N-B^ (lg) 
K su P A-B = V A-B K N~B^ (lg) 

where N-B refers to a nucleon nucleus collision. 

We know, without any calculation, that the Dual Parton Model for nucleus- 
nucleus collisions differs from the restricted superposition model. There are at 
least three reasons for these differences: 

1. Properties of the Glauber cascade: The deviations of the Glauber model 
from the restricted superposition model were already discussed in the previ- 
ous subsection. In Table | we give for each collision considered the average 
numbers of collisions according to the Glauber cascade, u, v p and u t were al- 
ready explained above. In nucleon-nucleus collisions we have always v N ~ B 
= vf~ B and v^~ B = 1. In A-B collisions we have u A ~ B , v A ~ B and v A ~ B 
all different from one. In the A-B collisions, a projectile nucleon interacts, 
in average, N A ~ B = v A ~ B jv A ~ B times. As seen from Table |3|, N A ~ B is 
always different from v^~ B at the same energy, N A ~ B being smaller than 
the corresponding vf~ B . The average nucleon-nucleus collision in an A-B 
collision has different properties compared to an average nucleon-nucleus 
collision. This is valid for any model constructed on the basis of the Glauber 
cascade. 

2. We turn to the DPM with Glauber cascade, and let us for simplicity dis- 
cuss the chain structure disregarding the multiple soft chains and multiple 
minijets in each nucleon-nucleon interaction, which are in addition present 
in each elementary Glauber collision. We construct in nucleon-nucleus col- 
lisions one pair of valence-valence chains and vf~ B - 1 pairs of sea-valence 
chains. In the restricted superposition model for the A-B collision we ob- 
tain v A ~ B pairs of valence-valence chains and v A ~ B {u^~ B — 1) pairs of 
sea-valence chains. 

In A-B collisions we have 

U A-B } V A-B 

and u A B all different from one and 
in Fe-Air collisions we have usually v A ~ B > u A ~ B > v A ~ B . In this situation 
we might form v A ~ B valence-valence chain pairs, v A ~ B — v A ~ B valence- 
sea chain pairs and v A ~ B — v A ~ B sea-sea chain pairs (this is the simplest 
possibility, in actual Monte Carlo models like DPMJET-II the chain ends 
available are connected in a random way). We know already that only 
the number of valence-valence chain pairs in the restricted superposition 
model and in the correct model correspond to each other, the number of 
sea-valence chain pairs will be different. Furthermore, there are never any 
sea sea chains in the restricted superposition model. 
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3. The single diffractive cross sections in hadron-nucleus collisions are well 



studied [^3], [34[] and included in the model. The cross sections for single 
diffractive nucleus-nucleus collisions are straightforward to calculate within 
the DPM, and they will probably have a size similar to that obtained in 
hadron-nucleus collisions. The O-air or Fe-air total cross sections are three 
to five times larger than the p-air total cross section. Therefore, the fraction 
of single diffractive nucleus-nucleus collisions will come out to be a few times 
smaller than in hadron-nucleus collisions. For these reasons, and because 
of the lack of experimental data on single diffractive nucleus-nucleus cross 
sections, this component is not included at present in the DPMJET-II event 
generator (however, the inclusion of this part would probably not change 
the situation). Now, in the case of diffractive p-air collisions, where the 
projectile proton leaves the collision only slightly deflected and in which 
hadrons are produced only in the nucleus fragmentation region, nearly no 
contribution to the moments of pions and kaons is achieved. Therefore, also 
this effect is expected to lead to larger moments of charged pions and kaons 
in the correct DPM with respect to the restricted superposition model. 

Looking at Tables p] and |2| we find our expectations fulfilled. The spectrum 
weighted moments of pions and kaons according to DPMJET-II are 10 to 20 % 
larger than the ones according to the restricted superposition model. Also the 
energy fractions carried by pions and kaons are for He-air and O-air about 5 to 
10 % larger than in the restricted superposition model. In Fe-air these energy 
fractions behave in a more complicated way. In the Fe-air reaction, it would be 
better to define the restricted superposition model in the Fe rest frame with the 
air nuclei as projectiles. At the same time, as expected, the energy fractions 
carried by the participating nucleons are for Fe-air and O-air in DPMJET-II 
always smaller than in the restricted superposition model, i.e. more energy is 
used for particle creation. Here the numbers obtained in He-air are more difficult 
to understand. We have to remember that all numbers in the Tables are the result 
of Monte Carlo calculations with typically 5% statistical error. 

In summary, we may conclude that for the model without FZIC we find an 
agreement with the restricted superposition model within approximately 10 to 
20 %. 



4.4 Calculation results, the full DPMJET-II model 

In Table |3| and |] we present spectrum weighted moments and energy fractions as 
calculated in the full DPMJET-II model. We consider again p-air, He-air, O-air 
and Fe-air collisions at laboratory energies per nucleon of 50 GeV, 500 TeV, and 
5 ■ 10 6 TeV. We also give the number of Glauber collisions v between v v projectile 
and v t target nucleons in minimum bias collisions in Table [3| 

In Tables |3] and |] we have a further comparison, at two energies, of the spec- 
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trum weighted moments of pions and kaons and of the energy fractions of charged 
pions, charged kaons, baryons minus antibaryons and of nuclear fragments and 
residual nuclei, with the moments obtained from the restricted superposition 
model. Here we have to use slightly different expressions for the energy fractions 
in the restricted superposition model. The reason is that all spectators in the 
full DPMJET-II appear among the final state particles, either as evaporation 
protons or neutrons or as nuclear fragments and residual nuclei. We are not able 
to define nuclear fragments and residual nuclei in the restricted superposition 
model. In this respect the features of DPMJET-II are again different from those 
of the restricted superposition model. We use the expressions for moments and 
energy fractions for charged pions and kaons in the same form as given in the 
last subsection. For the energy fraction of b — b we use: 

K1 A ~ B = <~ B K-f + {£* - <~ B ) ■ i- (20) 

Since we are not able to define an energy fraction K n _f, for nuclear fragments, we 
should compare the K^ A ~~ B in Table 4 with the sum of and K^ B . 

The differences between the spectrum weighted moments of the pions in the 
restricted model and in the full model in Tables [I] and ^ are easy to understand 
looking at the previously presented Feynman-x^ distributions in the projectile 
fragmentation regions. Comparing the results of DPMJET-II with full FZIC 
with the restricted superposition model in Tables ^| and |4] we find: 

1. For He-air collisions, and presumably for all light nuclei, the full model 
agrees practically to the model without formation zone cascade. Therefore, 
we also find only a 20 % disagreement with the restricted superposition 
model. 

2. We find that, for O-air and Fe-air, differences in the spectrum weighted 
moments of pions and kaons are up to 50 - 70% with respect to the restricted 
superposition model. We understand the reasons for this large disagreement 
looking at the changes in the Feynman xp distributions due to the FZIC as 
discussed in the last Section. The same trend is seen in the energy fractions 
of pions and kaons, but there the differences are only around 10%. 

3. The energy fractions of nuclear fragments K n .f. exist only in the full model. 
We find however the sum K b _i + K n j of the full model remarkably close 
to the energy fractions KT% of the restricted superposition model. 

5 Conclusions and summary 

We have compared Feynman-x^ distributions, spectrum weighted moments and 
energy fractions of pions, kaons and baryons minus antibaryons in two versions 
of the two component Dual Parton Model. One version is the full model with 
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FZIC of the produced hadrons with the spectators of the target and projectile 
nuclei, the second version is the DPM (closer to a restricted superposition model) 
without this formation zone cascade, nuclear evaporation and the formation of a 
residual nucleus. 

We have discussed the reasons for deviations expected between the DPM and 
restricted superposition models. Nevertheless, we find a reasonable, however not 
perfect, agreement of the moments from the restricted model for nucleus-nucleus 
collisions with a restricted superposition model. 

In the full model instead, the hadrons produced in the projectile fragmentation 
region by the formation zone cascade lead to significant deviations in the pion 
spectrum weighted moments from the restricted superposition model. For Fe-air 
collisions this difference becomes as large as 70%. This effect depends on the 
changes in the Feynman-a;^ distributions due to the cascade in the projectile 
nucleus. We have stressed above that there is no direct experimental evidence 
for this. Therefore we cannot be completely sure about the quantitative size of 
the predicted effect, but we are convinced that this effect exists. 

These differences are certainly significant and could also lead to differences 
when sampling the cosmic ray cascade with the full DPMJET-II model as com- 
pared to the sampling using the superposition model. Larger values for Z v and 
K w (also for 7r°, where the differences are similar, and kaons) mean that, on the 
average, the primary cosmic ray loses a larger fraction of its energy in the first in- 
teraction. This energy is spent for pion production. On the other hand, we expect 
that less energy is carried away by the leading particles, so that the development 
of the shower might be different from that resulting from the superposition model. 
Depending on the relevance of the first interaction, this should probably lead to 
showers with a larger muon content and with a smaller depth of the maximum 
development (thus with a smaller electromagnetic size at mountain or sea level). 
We might also expect that in the simulation of the full cosmic ray cascade initi- 
ated by primary nuclei, the full model could well show further deviations from the 
superposition model which are not simply contained in the differences of spec- 
trum weighted moments and energy fractions of pions and kaons. This is another 
reason to stress the importance of studying the full model in shower simulations. 

In summary, DPMJET-II including FZIC offers now the opportunity to per- 
form a true quantitative test of the validity (or failure) of the superposition 
approximation. A detailed study of the effects of the FZIC model on cosmic ray 
showers as a function of the primary energy and mass will be the object of a next 
paper. 
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Table 1: Spectrum weighted moments of pions and kaons according to 
DPMJET-II in the model without FZIC. In the last columns we compare with 
the restricted superposition model. E is the lab-energy per projectile nucleon. 



Collision 


E (TeV) 


Z n 


Z K 


u p 


v v Zl~ air 

P 7T 


ryp—air 
V V Zj K 


p-air 


500 


0.0640 


0.0092 


1 








5- 10 6 


0.0595 


0.0087 


1 






He-air 


500 


0.178 


0.025 


2.32 


0.148 


0.021 




5- 10 6 


0.180 


0.026 


2.60 


0.155 


0.023 


O-air 


500 


0.425 


0.059 


5.56 


0.356 


0.051 




5- 10 6 


0.462 


0.065 


6.49 


0.386 


0.057 


Fe-air 


500 


0.912 


0.124 


11.61 


0.743 


0.107 




5- 10 6 


1.020 


0.128 


14.46 


0.860 


0.126 
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Table 2: Energy fractions of pions , kaons and baryons minus antibaryons ac- 
cording to DPMJET-II in the model without FZIC In the last columns we com- 
pare with the restricted superposition model. E is the lab-energy per projectile 
nucleon. 



Collision 


E (TeV) 


K w 


K K 




v p 




v v K v K 


K p-_Air 


p-air 


500 


0.336 


0.048 


0.268 


1 










5- 10 6 


0.361 


0.052 


0.193 


1 








He-air 


500 


0.89 


0.13 


0.82 


2.32 


0.78 


0.11 


0.62 




5- 10 6 


1.01 


0.14 


0.76 


2.60 


0.94 


0.14 


0.50 


O-air 


500 


2.02 


0.28 


1.22 


5.56 


1.87 


0.27 


1.49 




5 • 10 6 


2.43 


0.32 


1.24 


6.49 


2.34 


0.34 


1.25 


Fe-air 


500 


4.05 


0.54 


2.37 


11.61 


3.90 


0.56 


3.11 




5- 10 6 


4.70 


0.60 


2.35 


14.46 


5.22 


0.75 


2.79 
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Table 3: Spectrum weighted moments of pions and kaons according to 
DPMJET-II in the model with full FZIC. We give also the values of is, v v and v t 
for the minimum bias Glauber cascade introduced in Section 2.1. E is the lab- 
energy per projectile nucleon. In the last columns we compare with the restricted 
superposition model. 



Collision 


E (TeV) 


V 


u p 






%k 


rjp-air 




p-air 


500 


2.05 


1. 


2.05 


0.0640 


0.0091 








5- 10 6 


2.73 


1. 


2.73 


0.0586 


0.0085 






He-air 


500 


4.44 


2.32 


3.41 


0.182 


0.026 


0.148 


0.021 




5- 10 6 


6.58 


2.60 


4.37 


0.175 


0.025 


0.152 


0.022 


O-air 


500 


10.34 


5.55 


5.32 


0.499 


0.065 


0.355 


0.051 




5- 10 6 


14.35 


6.22 


5.81 


0.522 


0.065 


0.364 


0.053 


Fe-air 


500 


21.51 


11.89 


6.87 


1.274 


0.137 


0.761 


0.101 




5- 10 6 


32.48 


13.86 


7.25 


1.287 


0.137 


0.812 


0.118 
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Table 4: Energy fractions of pions , kaons , baryons minus antibaryons and 
nuclear fragments according to DPMJET-II in the model with full FZIC. We 
give also the values of v v for the minimum bias Glauber cascade introduced in 
Section 2.1. E is the lab-energy per projectile nucleon. In the last columns we 
compare with the restricted superposition model. 



Collision 


E (TeV) 


u p 


K w 


K K 




K n .f. 


T^p-air 


v K v ~ air 


K sup 
b—b 


p-air 


500 
5- 10 6 


1. 
1. 


0.35 
0.37 


0.049 
0.053 


0.25 
0.18 










He-air 


500 


2.32 


0.90 


0.13 


1.58 


0.52 


0.81 


0.11 


2.26 




5- 10 6 


2.60 


0.99 


0.14 


1.37 


0.49 


0.95 


0.14 


1.87 


O-air 


500 


5.55 


2.22 


0.29 


4.35 


7.36 


1.93 


0.27 


11.84 




5- 10 6 


6.21 


2.56 


0.33 


4.11 


7.17 


2.27 


0.33 


10.91 


Fe-air 


500 


11.89 


4.81 


0.57 


18.01 


28.19 


4.14 


0.58 


47.08 




5- 10 6 


13.86 


5.43 


0.63 


16.03 


28.92 


5.06 


0.73 


44.63 
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Figure 1: Feynman-a;^ distribution of 7r + in p-air collisions at a lab-energy of 
500 TeV. Plot 1 (large full circles) is for the full model with FZIC in the target 
nucleus; plot 2 (small full circles) is for the model without FZIC. The target 
fragmentation region is at negative Feynman-Xi?, the projectile fragmentation is 
at positive xp- For comparison we give also the same distribution in p-p collisions 
(small empty circles). 
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Figure 2: Feynman-x^ distribution of 7r + in He-air collisions at a lab-energy 
of 500 TeV per projectile nucleon. Plot 1 (large full circles) is for the full model 
with FZIC in the target nucleus, and projectile nucleus; plot 2 (small full circles) 
is for the model without FZIC. The target fragmentation region is at negative 
Feynman-x^, the projectile fragmentation is at positive xp. For comparison we 
give also the same distribution in p-air collisions with FZIC (small empty circles). 
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Figure 3: Feynman-x^ distribution of ir + in O-air collisions at a lab-energy of 
500 TeV per projectile nucleon. Plot 1 (large full circles) is for the full model 
with FZIC in the target nucleus, and projectile nucleus; plot 2 (small full circles) 
is for the model without FZIC. The target fragmentation region is at negative 
Feynman-a;^, the projectile fragmentation is at positive xp. For comparison we 
give also the same distribution in p-air collisions with FZIC (small empty circles). 
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Figure 4: Feynman-Xi? distribution of 7r + in Fe-air collisions at a lab-energy 
of 500 TeV per projectile nucleon. Plot 1 (large full circles) is for the full model 
with FZIC in the target nucleus, and projectile nucleus; plot 2 (small full circles) 
is for the model without FZIC. The target fragmentation region is at negative 
Feynman-a;^, the projectile fragmentation is at positive xp. For comparison we 
give also the same distribution in p-air collisions with FZIC (small empty circles). 
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Figure 5: Glauber profile functions T(b) for Fe-nitrogen, He-nitrogen and p- 
nitrogen collisions at 500 TeV. 
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